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ABSTRACT 

A search for novae in M49 (NGC 4472) has been undertaken with the Hubble 
Space Telescope. A 55-day observing campaign in F555W (19 epochs) and F814W 
(five epochs) has led to the discovery of nine novae. We find that M49 may be 
under-abundant in slow, faint novae relative to the Milky Way and M31. Instead, 
the decline rates of the M49 novae are remarkably similar to those of novae in 
the LMC This fact argues against a simple classification of novae in "bulge" 
and "disk" sub-classes. We examine the Maximum- Magnitude versus Rate of 
Decline (MMRD) relation for novae in M49, finding only marginal agreement 
with the Galactic and M31 MMRD relations. A recalibration of the Buscombe- 
de Vaucouleurs relation gives an absolute magnitude 15 days past maximum of 
My_i5 = —6.36 ± 0.19, which is substantially brighter than previous calibrations 
based on Galactic novae. Monte Carlo simulations yield a global nova rate for 
M49 of rj = IOO^Iq year~^ and a luminosity-specific nova rate in the range uk = 
1.7 — 2.5 year~^ 1O~^°Lk0. These rates are far lower than those predicted by 
current models of nova production in elliptical galaxies and may point to a relative 
scarify of novae progenitors, or an increased recurrence timescale, in early-type 
environments. 

Subject headings: galaxies: individual (M49, NGC 4472) — novae, cataclysmic 
variables — galaxies: star clusters — stars: distances 



1. Introduction 

As close binaries in which material is accreted onto the surface of a white dwarf, novae 
form the cataclysmic variable (cv) subclass of variable stars. With amplitudes of 10-20 
magnitudes, they reach maximum magnitudes of —6.5 < My ^ —10 soon after the onset of 
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thermonuclear runaway burning. These high luminosities — coupled with their occurrence 
in galaxies of all morphological types — suggests that novae have the potential to be useful 
distance indicators, provided that some aspect of their behavior near maximum light can be 
used as a standard candle. 

While the earliest observations of extragalactic novae were reported by Ritchey (1917) 
and Shaplcy (1917), it was Hubble's exhaustive study of M31 that probably constituted the 
first dedicated search for novae in an external galaxy (Hubble 1929). The full potential of 
novae as distance indicators became apparent following the discovery by Zwicky (1936) that 
their peak brightness correlates with their rate of decline, in the sense that bright novae 
fade more rapidly than their faint counterparts. Although the nature and calibration of this 
"Maximum Magnitude versus Rate of Dechne" (MMRD) relation has been investigated on 
many subsequent occasions {e.g., McLaughhn 1945; Arp 1956; Cohen 1985; Capaccioh et al. 
1990; Livio 1992; Downes & Duerbeck 2000), these investigations have usually rehed on 
observations for a handful novae in the Galaxy or in M31, which, by virtue of its proximity 
and high luminosity, has remained the premier target for extragalactic novae surveys {e.g., 
Arp 1956; Rosino 1973; Rosino et al 1989; Shatter & Irby 2001). 

Despite some notable exceptions {e.g., Graham 1979; Pritchet & van den Bergh 1985; 
1987), studies of extragalactic novae have remained largely serendipitous in nature {e.g., 
Ferrarese et al. 1996). While a few heroic attempts to detect and study novae in Virgo 
and Fornax ellipticals using ground-based telescopes (Pritchet & van den Bergh 1985; 1987; 
Shafter, Ciardullo & Pritchet 2000; Delia Valle & Gilmozzi 2002) yielded light curves of 
varying quality for a few novae in a handful of galaxies, conclusions regarding the universality 
of the MMRD relation and its potential as a distance indicator rely almost entirely on 
observations of novae in the Galaxy, M31 and LMG {e.g., Delia Valle & Livio 1995). Accurate 
data for a sample of novae belonging to an elliptical galaxy with a well known distance would 
be invaluable in this regard, particularly since novae provide one of the few direct probes of 
compact binaries in such environments. 

In this paper, we report the results of the first dedicated search for extragalactic novae 
with the Hubble Space Telescope (HST). HST is an ideal instrument for such a survey thanks 
to its high spatial resolution, its ability to reach faint magnitudes in relatively short expo- 
sures, and the opportunity to schedule observations according to a pre-defined, optimized 
sequence. Our target, M49 (NGG 4472), is an obvious choice for several reasons. Not only 
it is the first ranked member of the Virgo Cluster, it is also the optically brightest galaxy 
in the local supercluster. Moreover, a variety of distance estimates are available for both 
Virgo, and for M49 itself. In fact, M49 offers the opportunity to compare, and perhaps even 
cahbrate, the nova MMRD relation directly against other Population II distance indicators 
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such as surface brightness fluctuations and globular clusters. 

2. Observing Strategy 

The if^r observations of M49 began on April 09, 2001 using the Wide Field and Plane- 
tary Camera 2 (WFPC2) on board the Hubble Space Telescope {HST). A total of 38 F555W 
(~ Johnson V) images, divided among 19 epochs, were obtained within a 55 day period, 
with the sequence ending on June 03, 2001. For epochs 14 to 18 of the sequence, F814W 
(~ Johnson J) images were taken immediately following the F555W observations. The rel- 
evant parameters of the observations are summarized in Table 1. All images were obtained 
at the same telescope roll angle. To facilitate removal of cosmic rays, exposures taken at 
each epoch shared the same pointing; however, the telescope was dithered by 0"6 between 
different epochs, to allow for accurate removal of CCD anomalies (e.g. hot or dead pixels). 
To maximize the number of novae outbursts, the center of M49 was placed close to the 
intersection of the four WFPC2 chips (Figure 1). 

The temporal sequence and duration of the exposures were chosen with two main consid- 
erations in mind, namely: the need to produce well-sampled novae light curves (particularly 
near maximum light) and the need to follow their decline for at least two magnitudes there- 
after. Because novae are not periodic (at least not on short timescales), there is always an 
obvious gain in lengthening the duration of monitoring programs. In practice, the length of 
our program was imposed by the pragmatic requirement of monitoring the field without a 
change in position angle, which HST can do for a maximum of about two months. This is 
longer than the decline rate for most novae (Capaccioli et al. 1989; Capaccioli et al. 1990; 
Downes & Duerbeck 2000), and therefore adequate for our purposes. 

Prior to outburst, nova progenitors are obviously much too faint to be detected at the 
distance of M49. Following outburst, they reach maximum light rather quickly (typically 
within a few days) and then decline over much longer periods: from several days to a few 
months. A two day interval between visits was deemed sufficient to sample the light curve 
adequately both before and immediately after maximum. In the interest of keeping the total 
time request for the program within reasonable limits, M49 was originally scheduled to be 
monitored every two days with the F555W filter for the first 14 epochs (26 days), and every 
five days thereafter for the final five epochs (25 additional days) with both the F555W and 
F814W filters. Although novae outbursts occurring during the second half of the sequence 
might not be optimally sampled, it was hoped that the time and magnitude at maximum, as 
well as the rate of decline, could be recovered with the aid of the color information (Ferrarese 
et al. 1996; see also §4.1) and by employing the well-sampled light curves of novae discovered 
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during the first part of the program as templates. 

The observing strategy discussed above was disturbed by a telescope safing event which 
occurred from 04/28/2001 to 04/30/2001. Because of this interruption, what was planned 
to be the eleventh epoch of our sequence was skipped and rescheduled at the end of the 
program, leaving an unfortunate four day gap during the finely sampled portion of the 
observing sequence. 

3. Data Analysis 

3.1. Data Reduction 

The WFPC2 consists of four separate 800x800 CCD detectors: three Wide Field Cam- 
era (WFC) chips, with a pixel size of 0.10 arcsec and a field of view of 1.3x1.3 arcmin per 
chip, and one high resolution Planetary Camera (PCI), with a pixel size of 0.046 arcsec 
and a field of view of 36x36 arcsec. The gain and readout noise are about 7 e~/DN and 
5 e~ respectively. Further details about the instrument can be found in the HST WFPC2 
Instrument Handbook (Biretta et al. 2001). All of the M49 observations were obtained with 
the telescope guiding in fine lock, which achieves a nominal pointing stability of about 3 
milliarcseconds. 

The reduction of the M49 frames followed the standard pipeline maintained by the Space 
Telescope Science Institute (STScI), and included correction of small A/D errors; subtrac- 
tion of a bias level for each chip; subtraction of a superbias frame; subtraction of a dark 
frame; correction for shutter shading effects and division by a fiat field. For each filter, two 
back-to-back exposures were obtained to aid in the removal of cosmic rays. After standard 
reduction, but prior to performing the photometric analysis, these exposures were combined, 
and cosmic rays fiagged and removed by comparing the difference in values between pairs 
of corresponding pixels to a local sigma calculated from the combined effects of Poisson 
statistics and local noise. The final, reduced, cosmic ray cleaned set of images consists of 
19 F555W and five F814W frames. As a final step prior to photometric reduction, the vi- 
gnetted edges of each chip were blocked and the geometric distortion of the WFPC2 optics 
was corrected using a pixel area map, following Stetson et al. (1998). 
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3.2. Photometric Analysis 

Photometric analysis of the data was performed using a variant of DoPHOT (Schechter 
et ah 1993), developed specifically to handle the peculiarities of the HST/WFPC2 Point 
Spread Function (PSF) (Saha et al. 1994). As is standard procedure (e.g. Saha et al. 1996), 
DoPHOT was first run on deep F555W and F814W "template" frames made by coadding all 
38 F555W and 10 F814W images, respectively. To prevent DoPHOT from triggering on an 
unreasonable number of noise spikes, a minimum S/N = 4 was required for object detection. 
The resulting master star list was then used as a position reference when DoPHOT was 
run on each individual frame. Using this approach, the completeness limit of the resulting 
photometry reaches fainter magnitudes than it would if DoPHOT were to be run directly 
on each individual epoch without the aid of a master star list. Furthermore, the results are 
more robust against residual cosmic rays which might affect individual frames. 

The major challenge posed to the photometry comes from the fact that the vast majority 
of the point-like objects in the WFPC2 field are globular clusters belonging to M49's rich 
system. At the distance of Virgo, globular clusters are marginally resolved, so that their PSF 
differs from the stellar PSF (and indeed varies with the size of the cluster). Because DoPHOT 
constructs its PSF directly from the brightest and most isolated objects in the image, in the 
case of M49 such a PSF might not adequately represent truly unresolved objects, such as 
novae. If this were the case, then the nova photometry could be compromised since the 
photometric corrections needed to convert the fitted DoPHOT magnitudes to conventional 
magnitudes (see below) would be appropriate for the globular clusters only. 

Fortunately, this appears not to be the case. For a few of the epochs, the fitted DoPHOT 
magnitudes (which are simply related to the height of the fitted PSF) were compared to those 
measured from test runs in which DoPHOT was forced to use a fixed PSF, constructed from 
HST/WFFC2 observations of the globular clusters Pal 4 and NGC 2419 (Stetson et al. 
1998). No significant offsets or scale errors were measured for either filter, indicating that 
the photometry is not affected by the small PSF broadening driven by the globular clusters. 

DoPHOT fit magnitudes differ from 'conventional' magnitudes by an additive aperture 
correction, which can be calculated as the difference between fit magnitudes and aperture 
magnitudes, provided that the latter can be obtained for a sufficient number of bright, iso- 
lated stars in the field. Since such stars are not available for M49, a more robust approach is 
to use aperture corrections measured from unrelated F555W and F814W WFPC2 observa- 
tions of a densely populated field in the Leo I dwarf galaxy (Saha et al. 1996). These images 
contain a multitude of stars across the entire field, allowing to quantify the positional depen- 
dence of the aperture corrections due to the spatial variations of the PSF. Conventionally, 
aperture magnitudes are measured within a 9" x 9" aperture; although it does not contain 
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the total light from a star, this aperture is large enough that PSF changes across the field, 
changes in focus over time, and small changes in jitter from one exposure to another do not 
affect the fraction of the total light that lies within this aperture. Finally, the mgxg magni- 
tudes obtained in this way were converted to the 'ground system' magnitudes F555W and 
F814W as defined in Holtzman et al. (1995) using the zero points derived from observations 
of ui Cen (Hill et al. 1998). In practice, this transformation is applied only to the results 
from the deep, template frames. Once the objects from the individual epochs are matched 
to their counterparts on the template, the offset for the individual epoch is evaluated from 
the ensemble average of magnitude differences (object by object) between that particular 
epoch and the template frame. 

Following Holtzman et al. (1995), F555W and F814W magnitudes can be converted to V 
and / by applying a correction which depends on the color of the object under consideration. 
While color information is not always available for the novae, this correction is always smaller 
than 0.02 mag in V and 0.03 mag in I over the color interval 0.2 < {V — I) < 0.7 mag, a range 
that is typical for novae (see, e.g., van den Bergh & Younger 1987). Although systematic, 
this error is insignificant compared to the other sources of errors which enter our analysis, 
and will be neglected for the remainder of this paper. 



Variability was searched for in the F555W frames using two independent methods. The 
first method applies a simple test to the DoPHOT photometry for each object. The 
second method is similar to the image subtraction technique described by Alard & Lupton 
(1998), and makes no use at all of the DoPHOT photometry. 

For any object detected by DoPHOT in at least two of the F555W frames, a reduced 
Xr was calculated as 



where rrii and tTj are the magnitude and rms error of a particular object as measured in 
the i—th epoch, m is the magnitude of the same object averaged over all epochs, and n is 
the number of epochs in which the object is detected. An object was fiagged as variable if 
> 3; this is effectively equivalent to selecting objects whose magnitude fiuctuates around 
the mean by 1.7a on average. 



3.3. 



Variable Star Seeirch 
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Most of the objects flagged are, in fact, not intrinsically variable; rather, the large 
is triggered by various anomalies in the images (e.g. residual cosmic ray events). A 
visual inspection of the light curve for each putative variable was sufficient to identify and 
reject these cases. The remaining candidates were visually inspected by blinking all of the 
individual frames against each other. The final list consists of nine bona-fide novae, all in the 
WFC chips (§4). Aperture photometry performed for each nova on all epochs in which the 
nova was detected provided an additional confirmation as to the rehability of the DoPHOT 
photometry. 

The image subtraction technique can be considered as a "visual" application of the 
method described above. First, all epochs are shifted to a common reference frame. A 
standard deviation frame is then created as: 

where fi{x,y) is the number of counts detected at pixel {x,y) in the i — th epoch F555W 
frame, and f{x,y) is the mean fiux at the same pixel, averaged over all epochs. To avoid 
hot pixels and/or single cosmic ray hits that inflate the standard deviation, at each pixel, 
the epoch with the largest value of fi{x,y) is excluded from the summation in equation 2. 
All objects detected at the 2 a level relative to the mean local background in the standard 
deviation frame were then visually inspected by blinking all of the original frames against 
each other. This procedure recovered all of the nine novae flagged by the test, but did 
not produce any additional detections. 

The location of the nine novae found in M49 is shown superimposed to gray scale images 
of each chip in Figures 2-4, and on an isophotal contour of M49 in Figure 5. Zoomed-in 
snapshots of the novae around maximum light are shown in Figures 6 and 7. The coordinates 
of the novae are given in Table 2, while Table 3 lists the magnitudes at each epoch. Light 
curves are shown in Figures 8-10. 

4. The Novae Light Curves 

Can we be sure that the nine variable objects are indeed novae belonging to M49? We 
immediately rule out the possibility that they are Galactic RR Lyrae or Cepheid variables on 
the basis of their distinctive light curves: none of the light curves is periodic, the amplitude 
of variations would be atypically large for either class of variables, and the variation time 
scale is too long for RR-Lyrae variables. Furthermore, RR-Lyrae variables or Cepheids 
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reaching apparent magnitudes of V > 22 would need to have distances > 250 kpc, placing 
them outside the Galaxy. Nor could the variables be lower main sequence flare stars in the 
Galactic halo since variability in such stars occurs rapidly, on timescales of minutes to hours 
(Hoffmeister, Richter & Wenzel 1984). The possibility that the variables are distant AGNs 
is likewise untenable: optical fluctuations of 1-3 magnitudes in QSOs and Seyfert galaxies 
typically happen on timescales of several months to years, rather than weeks (Peterson 2001). 
Furthermore, there is no evidence for a nonstellar appearance in any of the sources. 

Distant supcrnovae are also unlikely candidates. Over the redshift range 0.3 ^ z < 1.2, 
Type la supcrnovae will have magnitudes of 20.5 ^ V ^ 26.5 at peak brightness {e.g., 
Schmidt et al. 1998; Nobili et al. 2003), falling by ~ 0.7 magnitudes about 15 days later {e.g., 
Hamuy et al. 1996). Although such distant supcrnovae would resemble novae in M49, the 
number expected in our survey is far too low. The rate of faint. Type la supcrnovae (which, at 
these magnitudes, outnumber Type II supcrnovae by a wide margin) is ~ 160 deg~^ year~^ 
over the range 0.3 < 2; < 1.2 (Pain et al. 2002). Our WFPC2 fleld covers an area of 
5.66 square arcminutes and the survey duration is 55 days, giving a dismal ~ 0.04 distant 
supcrnovae expected in our survey. 

More hkely, the variables belong to M49 itself, as indicated by the fact that their spatial 
distribution (Figures 2-5) is not uniform across the field, but concentrates towards the center 
of the galaxy. Based on the light curves, variation amplitudes, and the fact that M49 is 
an E2 galaxy, we can again eliminate the possibility of long period variables or Ccphcids. 
Microlcnsing events of the sort described by Baltz & Silk (2000; 2001) are also unlikely. 
Scaling from the event rate observed in M87, we expect approximately one event per month 
for M49 (E. Baltz; private communication) with no stipulation on timescale. From Figure 2 
of Baltz & Silk (2000), we expect events of duration tpwHM ^ 5 days to comprise only a 
small fraction of the overall event rate. Moreover, microlcnsing light curves are expected 
to be symmetric around the peak time, which is not the case for any of the five objects for 
which our observations bracket the inferred maximum. We conclude that microlcnsing is not 
a significant source of contamination, leaving novae in M49 as the only realistic explanation 
for the observed fiares. 

Figure 11 shows a comparison between the cumulative distribution for the novae and 
the underlying galaxy light, measured in the band. As is apparent from Figures 2-5, 
the novae appear to be more centrally concentrated than the galaxy light: a KS test shows 
an 8% probability that they were drawn from the same distribution. Artificial star tests 
(§7) predict that only novae in the innermost ~ 5 arcsec of the galaxy will elude detection, 
and that beyond 17 arcsec even novae which have already decfined by two magnitude at the 
time of discovery will be recovered. Missing novae in the innermost region would of course 
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only exacerbate the contrast between the novae distribution and the underlying galaxy light. 
Pritchet & van den Bergh (1987) also found that the cumulative distribution of novae in 
M49 does not follow the underlying galaxy light, although in their case the novae seem to 
be under-represented in the outer (r > 40") parts of the galaxy. These somewhat puzzling 
results deserve further investigation in future surveys. 



4.1. Determination of Light Curve Parameters 

As will be discussed in detail in §6, the use of novae as distance indicators requires 
an accurate knowledge of their light curves, in particular, the maximum magnitude, Knax, 
the time of maximum light, tmaxj and the time taken to decline two magnitudes from peak 
brightness, ^2- Here we discuss the determination of these parameters for our sample of 
novae. 

Novae #1,2,3,8: Unfortunately, the light curves for all of these novae were already de- 
clining from maximum brightness at the outset of our program, making direct measurements 
of Knax and tinax impossiblc. Instead, we take Knax ^V{1) and tmax < ^(1)) where V{1) and 
t{l) refer to the magnitude and time of the first epoch. The time taken for the novae to 
decline in brightness by two magnitudes, t2, and the corresponding decline rates, = 2/^2 
can be bracketed more securely, since the observed light curves are very nearly linear for the 
first few epochs. Both are therefore estimated directly from weighted least-squares fits to 
the K-band light curves using epochs 1-5 (nova #1), 1-6 (nova #2), 1-3 (nova #3), and 
1-10 (nova #8). 

Nova #4: The nova is detected both before and after maximum. Although the entire 
light curve is only sampled at 5-day intervals, both V and / data are available post-maximum, 
allowing a precise determination of tmax- According to van den Bergh & Younger (1987), 
novae at maximum light have {B — V)q = 0.23±0.06, corresponding to a spectral type 
between A7 and FO. From Zombeck (1990), this translates to {V — /)o = 0.37±0.10, or 
{V — /)max = 0.40 ± 0.10 in the case of M49, adopting a total (foreground plus internal) 
reddening to the galaxy A(y ) = 0.07 and the relative extinction coefficients from Schlegel, 
Finkbeiner & Davis (1998). Linear extrapolations from F(14), 1^(15), /(14) and /(15) then 
give tmax — 2452039.04 ± 1.11. The imphed magnitudes at this time are T4iax = 23.33 ± 0.13 
and /max = 22.93±0.03, where the main contribution to the error comes from the uncertainty 
in the intrinsic color at maximum. This nova has not yet faded by two magnitudes by the 
end of our observing sequence, but a weighted linear least-squares fit to V^(14)-\^(19) — 
further constrained to pass through V^ax — gives t2 — 52. 00 ±4. 21 days. 
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Nova #5: As in the preceding case, the nova is detected during its rise to maximum, 
and V and I Ught curves are available during its post-maximum evolution. Extrapolating 
to {V - /)max = 0.40 ± 0.10 on the basis of ^(15), ^(16), /(15) and J(16) gives = 
2452041.42 ± 1.95, Kiax = 23.10 ± 0.52 and I^^^ = 22.70 ± 0.42. A least-squares fit to V^^, 
V{IA) and 1^(15) yields ^2 = 7.53 ± 1.14 days. 

Nova #6: Although this nova is observed during its rise to maximum, only limited color 
information is available during its post-maximum evolution. We assume that peak brightness 
occurred midway between epochs 8 and 9, so that tmax = 2452023.46 ± 0.50. Extrapolations 
based on V^(9), V^(IO) and V^(ll) then yield V^max = 23.33 ± 0.28. A weighted least-squares 
fit to V{'3t)-V{lh) gives = 19.30 ± 0.86 days. 

Nova #7: Two observations of this nova are available prior to maximum light: V^(9) and 
y(10). In an unfortunate example of Murphy's law, the period around maximum coincided 
with the telescope safing event discussed in §2 (the light curve of nova 9 also suffered from 
the safing). Extrapolating forward from 1/(9) and V{1Q), and backwards from V{11) and 
1/(12), gives imax = 2452027.33 ± 0.50 and Kiax = 22.72 ± 0.36. Based on this V^^, and 
interpolating between 1/(13) and 1/(14), we find ^2 = 11.68 ± 0.83 days. 

Nova #9: In this case, there is some ambiguity concerning the precise time of maximum. 
Given the short rise time for most novae, we think it most likely that maximum brightness 
occurred shortly after epoch 11. The other possibility — that maximum brightness occurred 
before epoch 11 and after epoch 10, when the nova was not detected — cannot be excluded 
given the fact that novae light curves do not always exhibit a monotonic decline after max- 
imum (see, e.g., nova #6). Working under the assumption that peak brightness occurred 
midway between epochs 11 and 12, at tmax = 2452031.45 ± 0.50, extrapolating from 1/(12) 
and 1/(13) gives Knax = 23.74 ±0.22. This nova never fades by two magnitudes from Vmax so 
we take t2 = 37.16±2.86 days based on a linear least-squares fit to 1/(12)~1/(19), constrained 
to pass through Knax- 

Light curve parameters for all nine novae are summarized in Table 6. From left to 
right, the columns of this table record the identification number, the time of maximum light, 
the l/-magnitude at maximum, the absolute V magnitude for (m — M) = 31.06±0.06 and 
A(V)g = 0.07 (see §6.1), the time taken for the light curve to decline by two magnitudes, 
and the corresponding decline rate. 
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5. The Distribution of Novae Decline Rates 

It has been suggested on a number of occasions {e.g., Duerbeck 1990; Delia Valle & 
Duerbeck 1993; Delia Valle 2002) that there exist two populations of novae: bright, rapidly 
declining "disk" novae and faint, slowly declining "bulge" novae. Since M49 is an elliptical 
galaxy consisting mainly of old and intermediate- age stars {e.g., Trager et al. 2000; Cohen 
et al. 2003), one would expect our sample of novae to consist entirely of faint objects with 
long decline times. The cumulative fraction of novae as a function of decline rate, — 2/t2 
is shown in Figure 12, compared to the cumulative rates (again measured as Vd = 2/t2) 
observed for M31 (Capaccioli et al. 1989), the LMC (Capaccioh et al. 1990) and the Milky 
Way (Downes & Duerbeck 2000). The M49 sample seems to behave exactly the opposite as 
expected: i.e., faint objects with long decline rates are under-abundant. This is unlikely due 
to an observational bias since our photometry is essentially complete down toV = 25.5 mag. 
Even the faintest novae are expected to reach 23.8 < V < 24.5 mag at maximum (see §6.2) 
and decline by one magnitude or more during the length of our observing period: if present, 
they should have been detected. A Kolmogorov-Smirnov test yields significance levels of 
21% and 49% for the null hypothesis that the M49 novae sample and those of the Milky Way 
and M31, respectively, are drawn from the same population. However, the corresponding 
significance level for the M49 and LMC samples is 99%. 

It is worth noting that some of the decline rates for LMC and Milky Way novae, and 
most of those for novae in M31, are based on photographic magnitudes. Since novae decline 
more slowly in B than V, the decline rates for the LMC, Milky Way and M31 novae should be 
increased somewhat relative to those in M49, which are measured in the V band. According 
to Van den Bergh and Younger (1987), logt2(V^) = 0.953(±0.013) logt2{B). Correcting the 
Milky Way and LMC samples produces no difference as far as the KS test is concerned. For 
the M31 sample, it is not always clear which novae have photographic magnitudes and which 
have V— band magnitudes. In the extreme case in which the decline rates for all novae are 
corrected, the disagreement between the M49 and M31 samples becomes less severe, but a 
KS test still returns a significance level of only 84% that the two samples are drawn from 
the same population. 

The LMC novae are usually regarded as prototypical examples of young, disk novae. 
We conclude that our survey provides no clear support for a simple classification of novae 
into a disk and bulge populations, although a larger sample of novae in additional early type 
galaxies is needed for a clear resolution of the issue. 
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6. Novae as Distance Indicators 

The main goal of our program is to assess the usefulness of novae as standard candles 
and, if possible, to provide a calibration of the various relations which could establish novae 
as reliable distance indicators {e.g., Livio 1992). Because such relations are, at present, 
calibrated exclusively with novae in two spiral galaxies (i.e., M31 and the Galaxy), their 
applicability to early type galaxies remains untested. M49 was chosen for this program 
not only because of its brightness, but also because of the availabihty of secure distance 
estimates. Here, we focus on two Population II distance indicators — surface brightness 
fluctuations (SBF) and the globular clusters luminosity function (GCLF). 

6.1. Distance Calibration via SBF and Globular Clusters 

Tonry et al. (2001) report an SBF distance modulus of (m-M)o = 31.06±0.10 for MA9'^. 
Since the best available Maximum Magnitude versus Rate of Dechne (MMRD) relation for 
novae is based almost entirely on observations of M31, it is also of interest to know the 
distance to M31 on the same SBF scale, to facilitate the comparison of novae in the two 
galaxies. According to Tonry et al. (2001), the SBF distance of M31 is (m — M)q = 
24.40±0.08, giving a difference in distance moduli between the two galaxies of A(m — M)o 
= 6.66±0.13 mag. 

The GCLF has been claimed to be a first class distance indicator {e.g., Jacoby et al. 
1992) although concerns have also been raised (Ferrarese et al. 2000). Our co-added images 
of M49 are the deepest ever obtained for this galaxy; outside of the central 20", our deep 
images reach ~ 25 with better than 90% completeness. The GCLF, (^{V), can therefore be 
derived directly from our data, providing an independent check of the adopted SBF distance. 

To do so, the co-added, cosmic ray cleaned, deep F555W and F814W images (§3.1) 
were background-subtracted via multi-resolution wavelet filtering using the MR/1 package 
(Starck, Bijaoui & Murtagh 1998). Object detection and photometry on the background 
subtracted frames was performed with SExtractor (Bertin & Arnouts 1996), using a detection 
threshold of three connected pixels above 3.5a. Object catalogs for the separate F555W and 
F814W frames were then matched with a matching radius of two pixels, and calibrated 



^The SBF calibration adopted by Tonry et al. is based on Cepheid distances to nearby galaxies which 
use the same HST/WFPC2 photometric calibration adopted for the novae in this paper. Although the 
photometric scale of the WFPC2 has since been refined (Stetson 1998), the novae photometry and SBF 
distance are mutually consistent. 
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V and / magnitudes were obtained following Holtzman et al (1995, see also §3.2). The 
matched catalog was then trimmed to exclude objects with y < 25 or colors outside the 
range 0.7 < (V^ - /) < 1.45. 

Our final catalog consists of 389 globular cluster candidates. For comparison, Larsen 
et al. (2001) have recently carried out a study of the M49 globular cluster system based on 
three partially overlapping fields. Using the same magnitude and color selection criteria as 
described above, we find their catalog to contain a total of 661 objects; the larger sample 
is a consequence of the greater areal coverage of their survey. A comparison of the M49 
GCLFs derived in this paper and Larsen et al. (2001) is given in Figure 13 (note that our 
luminosity function has been scaled upwards by the ratio of the sample sizes, 661/389 ~ 
1.7). There is good agreement between the two luminosity functions. Following the usual 
procedure of parameterizing (f>{V) as a Gaussian, we find best-fit values for the turnover 
and dispersion to be = 23.87 ± 0.06 and a = 1.43 ± 0.09. With the calibration of 
Harris (2001), this turnover corresponds to a distance modulus of {m — M) = 31.13±0.09. 
This best-fit Gaussian is shown by the solid curve in Figure 13. If we instead choose to fix 
the location of the turnover based upon the SBF distance modulus of Tonry et al. (2001) 
and the M^o calibration of Harris (2001), we find = 23.80 and a = 1.42 ± 0.07. The 
corresponding luminosity function is shown as the dashed curve in Figure 13. We conclude 
that the SBF and globular cluster luminosity function methods yield highly consistent results 
for the distance of M49. 

In what follows, we adopt the SBF M49 distance modulus of (m — M) = 31.06 ± 0.10 
mag, corresponding to a distance of 16.3±0.7 Mpc. 

6.2. The Maximum Magnitude versus Rate of Decline Relation 

Although a relation between the magnitude reached by novae at maximum and their 
rate of decline was first proposed in 1936 by Zwicky, a theoretical and phenomenological 
description of this MMRD relation has proven remarkably elusive. From an observational 
standpoint, the main obstacle remains the small number of galaxies for which large samples 
of novae with high-quahty light curves have been collected. The hst includes only three 
objects: M31, with Ki 55 novae (Capaccioh et al. 1989), the Galaxy, with a few dozen 
objects (Downes & Duerbeck 2000), and the LMC, with about 10 novae (Capaccioli et al. 
1990). In many cases, analysis of the light curves is further complicated by the fact that 
observations are available in non-standard, photographic bandpasses. 

These concerns notwithstanding, it is known that the Galactic MMRD relation can be 
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fit reasonably well with a power law (Cohen 1985). Based on new distances derived from 
expansion parallaxes for a sample of Galactic novae, Downes & Duerbeck (2000) find: 

Mv,max = (-11.32 ± 0.44) + (2.55 ± 0.32)log(i2) (3) 

where t2 is the time it takes for the nova to decline by two magnitudes from peak bright- 
ness, My^max- Both the M31 sample, and the combined M31 and LMC samples, are better 
described in terms of a "stretched" S-shaped curve, whose most recent calibration is given 
as (Delia Valle & Livio 1995): 

1.32-logi2 . X 

Mv,max = -7.92 - O.Slarctan (4) 

Fitting a similar function to the sample of Galactic novae produces a zero point consis- 
tent with the one found for M31, but a significantly higher contrast between faint and bright 
novae (Downes & Duerbeck 2000), although it should be pointed out that the quahty of the 
fit is not significantly improved over that provided by the single power law given in equation 
3: 

1.32-logt2 . . 

Mv,max = -8.02 - 1.23arctan q^S ' 

It is important to note that there is no theoretical explanation for a relation of the kind 
shown in equations 4 and 5. In fact, Downes & Duerbeck (2000) argue that a better, and 
physically motivated, characterization of the Galactic MMRD relation might be expressed 
in terms of a broken power-law, in which novae are divided in two distinct subgroups on 
the basis of the shape of the light curve and the detailed physics of the outburst (Duerbeck 
1981). 

Figures 14 and 15 compare the location of the nine M49 novae in the MMRD plane 
(circled points) with those of M31 and the Milky Way, respectively. Absolute, dereddened 
magnitudes and decline rates for the Milky Way novae are listed by Downes & Duerbeck 
(2000, Table 5), and can be easily scaled for comparison with the M49 sample by adopting 
(to - M)o = 31.06 ± 0.10 mag and A{V)g = 0.07 for M49 (from Tonry et al. 2001, using 
the relative extinction coefficients of Schlegel, Finkbeiner & Davis 1998). Dechne rates and 
apparent, reddened, photographic maximum magnitudes for the M31 novae are from Table 
IV of Capaccioli et al. (1989). The photographic magnitudes can be transformed to the 
l^— band according to the relation V ~ rripg — 0.06, but scaling the M31 novae samples to 
the distance and reddening of M49 proves to be a non-trivial task. Although the distance 
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to M31 is well known (e.g. Ferrarese et al. 2000, we adopt here (m - M)o = 24.40±0.08 
mag from Tonry et al. 2001), the total extinction to the galaxy is not. The RC3 reports 
values of A{B)i = 0.67 mag for the 5-band extinction internal to M31, based on the galaxy 
inclination and Hubble type. The foreground (Galactic) reddening is generally taken as 
E{B — V") = 0.08 mag, corresponding to A{B) = 0.35 mag (see Bianchi et al. 1996 for a 
comprehensive discussion). Given these estimates, the total K-band extinction to M31 is 
therefore A{V) ~ 0.77 mag, where we have adopted A{B) / A{V) ~ 1.321 following Schlegel, 
Finkbeiner & Davis (1998). This value can be compared to the one derived by Bianchi 
et al. (1996) using HST UV observations of eight blue supergiants. Their study points 
to a total reddening to M31 in the range E{B — V) — 0.13 — 0.20 mag, which translates 
to a band total extinction in the range A(y) — 0.43 — 0.66. For consistency with the 
absorption adopted for other galaxies with measured novae rates (to be discussed in §7), in 
what follows we adopt A{V) ~ 0.77 mag for M31, with a 20% uncertainty. 

Filled and open squares in Figure 14 represent novae judged by Capaccioli et al. to 
have light curves of good and fair quality, respectively. Absolute magnitudes are shown on 
the right vertical axis. The dashed curves in Figures 14 and 15 represent the best fit to the 
M31 sample (from Capaccioli et al. 1989), and the MW sample (Downed & Duerbeck, 2000; 
equation 5) respectively. Residuals between these relations and the M49 novae are shown 
in the lower panels. The sohd line in Figures 15 is the best power law fit to the Galactic 
sample (equation 3). 

The M49 sample is formally consistent with the Milky Way sample, while the agreement 
with the M31 data is less secure. Restricting ourselves to the five M49 novae with measured 
Viaax and Pd: the mean offsets between M49 and the best fitting relations to the M31 and 
Milky Way samples shown by the dashed lines in Figures 14 and 15 are —0.50 ± 0.15 
(random) ± 0.20 (systematic) mag and —0.18 ± 0.15 (random) ± 0.10 (systematic) mag 
respectively. In the above estimate, the systematic error arises from the uncertainties in the 
distances and reddenings to M31 and M49, while the random component arises from the 
uncertainties in the magnitudes at maximum for the M49 novae. Although the M49 and 
M31 samples disagree at the 2(T level, larger uncertainties in the internal extinction to the 
galaxy than considered here are possible. The root-mean square scatters of the M49 relation 
around the M31 and MW MMRD relations are 0.74 mag and 0.85 mag respectively; in both 
cases, there appears to be a systematic trend for fast novae to be under-luminous, and slow 
novae to be over-luminous compared to expectations, although the small number of objects 
does not allow us to speculate any further. 

A formal fit to the M49 sample, using the same relations given in equations 4 and 5, 
allowing only the zero point to vary, and accounting for errors in the magnitudes, produces: 



-16- 



Mv,max = (-8.27 ± 0.59) - O.Slarctan^ —; rms = 0.82 mag (6) 

0.23 

and 

1.32 — loEfto / s 

Mv,max = (-8.56 ± 0.86) - 1.23arctan ; rms = 1.24 mag (7) 

respectively. The large scatter about these relations suggests that, regardless of the details 
of the calibration, caution is needed when deriving distances via the MMRD relation. 



6.3. The Buscombe and de Vaucouleurs Relation 

Buscombe & de Vaucouleurs (1955) first noted that the absolute magnitude 15 days after 
maximum light, M15, appears to be roughly the same for all novae. The existence of a time, or 
time interval, after maximum at which the absolute brightness of novae are approximatively 
constant follows from the MMRD relation {e.g., Shara 1981b) — because brighter novae 
decline faster than fainter ones, the light curves must intersect sometime after maximum. 
The calibration of M15 as a standard candle has been attempted on many occasions {e.g., 
Shara 1981a; van den Bergh & Younger 1987; Cohen 1985; Downes & Duerbeck 2000), 
yielding a rather wide range in results. 

Figure 16 shows light curves for the five novae in our sample {^ 4, 5, 6, 7 and 9) for 
which both Knax and tmax can be measured directly from their light curves (see §4.1). Each 
light curve has been shifted in time so that maximum hght occurs at t = 0. Averaging the 
apparent magnitude at i = 15 days gives V15 — 24.77 ± 0.19, corresponding to 

Mv,i5 = -6.36 ± 0.19 (random) ± 0.10 (systematic) (8) 

where the random uncertainty represent the error in the mean, and the systematic uncer- 
tainty reflects the assumed errors on M49's distance and extinction. The standard deviation 
about this mean value is cr = 0.43. 

This determination of My^w is strongly at odds with the value My^w — —5.23 ± 0.16 
found by van den Bergh & Younger (1987) from a sample of six Galactic novae with known 
expansion parallaxes. It is barely consistent with the value of Mv,i5 — —5.60 ± 0.45 found 
by Cohen (1985) from her study of 11 novae with spatially resolved expansion shells, but 
it is in reasonable agreement with the value of My 15 = —6.05 ± 0.44 found by Downes & 
Duerbeck (2000) using the same method, but for a larger number of Galactic novae. 
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The lower panel of Figure 16 shows the mean novae magnitudes and standard deviations 
as a function of time after maximum. It is interesting to note that a rather small standard 
deviation [a = 0.37 mag) is found at maximum. The mean magnitude at maximum for the 
five observed novae is My^max = —7.89 ± 0.17 (random) ±0.10 (systematic), although we 
strongly caution against using this estimate as a distance indicator. For instance, nova ^^2 
was pa 0.6 magnitudes brighter than this at the onset of our program, at which time it was 
already declining from maximum. 

Is it reasonable to use My^i^ as a distance indicator? We can test this hypothesis on 
novae 7^1, 2, 3 and 8, which were already declining at the beginning of our program. If we 
assume that these novae reach V = 24.77 ± 0.19 fifteen days after maximum, then using the 
available datapoints (Table 3), we can linearly extrapolate the light curve and infer My^^ax 
and tmax- The V magnitudes at maximum found in this way are 20.87 ± 0.46 for nova # 1, 
21.53 ± 0.20 for nova # 2, 17.86 ± 1.1 for nova # 3, and 23.08 ± 0.19 for nova # 4^ If these 
data are added to the maximum magnitude and rate of decline of the five novae for which 
both parameters are measured directly, and the MMRD relation given in equation (5) is fit 
to the data, a distance modulus to M49 of (m — M)o = 30.67 ± 0.71 is obtained (the error 
accounts only for the uncertainties in the M49 data, not in the MMRD calibration). This is 
certainly consistent with Tonry et al. (m — M)q = 31.06 ± 0.10, but hardly useful. 

We conclude that the large standard deviation shown for Myjs, coupled with the rather 
poor agreement between the different calibrations of Myi5 based on samples of Galactic 
novae, suggest that one should exercise considerable caution when estimating distances with 
this method. This perhaps should come as no surprise since, as noted by Jacoby et al. (1992), 
several "exceptional" novae in M31 {e.g., Arp 1, 2 and 3) offered clear counterexamples to 
the universality of the Buscombe & de Vaucouleurs relation. 



7. Nova Rates 

The nova rate is a fundamental property of any stellar population and represents a 
direct probe of the abundance of compact, mass-transfer binaries in the host galaxy. There 
have been some recent claims that the luminosity-specific nova rate varies along the Hubble 
sequence, with late-type galaxies being more efficient producers of novae than their early- 
type counterparts {e.g., Delia Valle et al. 1994; Yungelson, Livio & Tutukov 1997; Delia 
Valle 2002; c.f. Shatter et al. 2000). M49, as an E2 galaxy and the optically brightest 
member of the local supercluster, presents an opportunity to measure the nova rate in an 



^For this nova, the outburst is predicted to have occurred 1.8 ± 2.7 days after the first epoch. 
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extreme environment. 

We use a Monte Carlo approach to measure the global nova rate in M49. The difference 
in distance moduli between M31 and M49 is taken to be A(m — M) — 6.66, as discussed 
in §6.1. For an assumed global nova rate, rj, we randomly select novae having maximum 
magnitudes and rates of dechne given by the M31 relation of Capaccioh et al. (1989), 
displaced to the distance of M49. Light curves are then simulated by calculating l^-band 
apparent magnitudes at the dates of the actual observations, assuming a linear decline from 
maximum light. 

For each simulated nova, we randomly assign a galactocentric position assuming that 
the local nova rate correlates linearly with stellar mass density, and adopting the M49 surface 
brightness model presented in Cote et al. (2003). Any novae falling outside the regions of 
our survey are discarded. Artificial star experiments were carried out to determine the level 
of photometric completeness as a function of radius and magnitude. The results of these 
experiments were used to determine, on a case-by-case basis, if the simulated nova would 
be detected in our survey at the assigned radius and given its instantaneous magnitude at 
each epoch. We consider a simulated nova to be "discovered" if it is detected at two or more 
epochs and shows a variation in brightness that is larger than that expected on the basis 
of its photometric errors. Wc carry out 1000 simulations at each t], recording the number 
of times the number of detected novae in the simulations matched the observed number of 
novae. This process is repeated for global nova rates varying between 77 = and 300 year~^ 
in steps of five. 

The results of this exercise are plotted in Figure 17. The histogram shows the number 
of matches found in the Monte Carlo simulations, plotted against global nova rate. The 
vertical lines outline the best estimate of the global nova rate, r] = 100^30 year^^, where the 
quoted uncertainties refer to 68% confidence limits determined directly from the simulations. 

How does this nova rate compare to those found in the literature? In a ~ 30 day 
campaign, Pritchet & van den Bergh (1987) visually identified eight novae in M49. Not 
accounting for incompleteness or detection biases, this yields a nova rate of 77 = 160 ± 
57 year~^ (for (m — M)q = 31.0), consistent with the value derived in this paper. The 
only other (undisturbed) giant elliptical galaxy for which a detailed search for novae has 
been undertaken is M87, the second brightest member of the Virgo cluster. On the basis of 
multi-epoch, ground-based Ha imaging, Shafter et al. (2000) measured a global nova rate of 
rj — 91 ± 34 year~^ for this galaxy. Scaling our nova rate for M49 downwards by the ratio 
of the K-hand luminosities of M49 and M87, we predict 77 ^ 70 ± 23 year~^ for M87, fully 
consistent with the nova rate of Shafter et al. (2000). On the other hand, our predicted 
nova rate for M87 is strongly at odds with the extreme rate of "200 to 300" year~^ found 
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for M87 by Shara & Zurek (2002)^. For instance, if we turn the argument around and scale 
the Shara & Zurek (2002) nova rate for M87 to the luminosity of M49, we would expect 
280 ^ 1] < 425 ycar~^. Our Monte Carlo simulations reveal that such high nova rates can 
be ruled out at better than 99.9% confidence. 

The complete sample of galaxies having measured nova rates encompasses objects with 
widely different morphological types, luminosities and star formation histories. Relevant 
information for these galaxies is summarized in Table 7. Columns (1-3) give the galaxy 
name, morphological type from the NASA Extragalactic Database (NED) and global nova 
rate. Sources for these nova rates arc given in the footnotes to the table. Note that there is 
some disagreement over the global nova rate in M33: Delia Valle et al. (1994) quote a value 
of 77 = 4.6 ± 0.9 year~^, while Sharov (1993) argues for an upper limit of 77 ~ 0.45 year~^. 

We follow the usual approach of normalizing the global nova rates by K-hand luminosity 
(which more closely traces the mass in evolved stars). We consider two estimates for the total, 
dereddened K-band magnitudes Mk,o- A recent 2MASS release (Jarrett et al. 2003) lists 
total (extrapolated), reddened, X-band magnitudes for all galaxies (Column 5 of Table 7), 
with the exception of the LMC and SMC. We corrected the 2MASS magnitudes for both 
internal and Galactic extinction, using the S-band absorption A{B)i and A{B)g from the 
RC3 and Schlegel, Finkbeiner & Davis (1998) respectively (columns 7 and 8 of Table 7), and 
A{K)/A{B) = 0.085 (Schlegel, Finkbeiner & Davis 1998). Mk,o can also be calculated by 
correcting the total 5-band magnitude (from the RC3, column 4 of Table 7) for extinction, 
and then applying a mean, dereddened (B — K)q color (column 6 of Table 7, with references 
given in the footnotes). To transform the dereddened Mk,o to a luminosity in solar units, we 
adopt the distance moduli given in column (9), and Mk,q — 3.33 mag (Cox 2000). Columns 
(10-11) give the luminosity-specific nova rates, i/k, in the two cases in which the galaxy K- 
band luminosity is calculated starting from the RC3 S-band magnitudes, or adopted from 
the 2MASS analysis, respectively. 

Before proceeding with an analysis of the behavior of among the full sample of 
galaxies, we pause to consider the implications of our measured nova rate for M49. Our best 
estimate of the luminosity-specific nova rate in M49, uk — 1.71±0.61 year~^ IO^^^Lk© {i^k — 
2.52 ± 0.91 year~^ 1O~^''Lk0 when using the 2MASS data), is far smaller than predicted by 
theoretical models of nova production. For instance, Matteucci et al. (2003) have computed 
galactic nova rates as a function of mass and star formation history. Using the elliptical 
galaxy template models, and adopting a luminous mass oi Ai ^ 8 x lO^^TVlo for M49 (Cote 
et al. 2003), we find the Matteucci et al. (2003) models to predict global nova rates in 



''As cited in Matteucci et al. (2003). 
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the range 780 ^ rj < 900 year ^, nearly an order of magnitude larger than the observed 
value oir] = lOOtH year . Given the size of the discrepancy — and since the theoretical 
models were calibrated to reproduce the global nova rate of 77 ~ 25 year"^ for the Milky Way 
— reconciling the models and observations for M49 may require global differences in the 
properties of the novae progenitors in the two galaxies. Obvious candidates would include 
a reduced binary fraction in M49 relative to the Milky Way, or longer recurrence timescales 
between outbursts. 

Figure 18 shows uk plotted as a function of K-hand luminosity for M49 (circled point) 
and each of the galaxies in Table 7 (dots). Note that M33 is plotted for both the high 
nova rate of Delia Valle et al. (1994) and the upper limit of Sharov (1993). Excluding this 
galaxy, the weighted mean nova rate is {uk) ~ 1-58 ± 0.16 year~^ 1O~^°Lk0 when AIk,o is 
derived from the RC3 S-band magnitudes, and (uk) — 2.41 ± 0.27 year~^ IO^^'^Lk© when 
the 2MASS data are used.^ 

The ~ 35% increase in the luminosity-specific nova rates when the 2MASS data are 
used is a consequence of the 0.2 mag systematic difference between the two estimates of 
Mk,o: with the 2MASS magnitudes being fainter. It is unclear which estimate of Mk,o is 
more reliable: compared to the 2MASS values, the RC3 S-band magnitudes are more robust 
against small errors in the sky estimates, but are very sensitive to extinction corrections, 
which can be very uncertain in the case of spiral galaxies. In either case, the luminosity 
specific nova rate in M49 appears to be in good agreement with the mean determined from 
the complete sample. Indeed, with the possible exception of the LMC, every galaxy for 
which a reliable measurement of uk exists appears to be consistent with a "universal" value 
of (uk) ~ 1.6 — 2.4 year"-*^ 10~^°Lko. These findings are fully consistent with those of Shaffer 
et al. (2000). 

This constancy of has implications for theoretical models of the galactic nova rates, 
and for possible dependences on morphological type and star formation history. Calculations 
by Yungelson, Livio & Tutukov (1997) predict a ~ 20-fold increase in the nova rate of a stellar 
population which has formed stars continuously over a Hubble time, compared to that of a 
stellar population which arose in a single burst 15 Gyr ago. While the precise star formation 
history of M49 is unknown, we note that this galaxy is likely to have a significant number 
of stars younger than 5-8 Gyr (Trager et al. 2000; Cohen, Blakeslee & Cote 2003). Thus, 
if we scale our luminosity-specific nova rate for M49 to late-type systems like the LMC and 
SMC, we would expect ^ 30 year~^ 1O~^°Lk0 for these galaxies. Microlensing surveys of 



^For comparison, the unweighted mean value are (vk) = 2.20 ± 0.24 year ^ 10 ^"Lk© and (i'k) = 
2.66 ± 0.30 year-i IO'^'^Lkq respectively. 
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these galaxies show that such extreme rates can be ruled out with high confidence, despite 
the large uncertainties in the measured values of 77 and vk- 

8. Conclusions 

We have presented the results of an if5'T'/WFPC2 program designed to discover novae 
in M49. Nine novae, five of which with fairly complete (i.e., covering both the pre- and post- 
maximum phases) and well-sampled light curves, were discovered in a 55-day campaign. 
These nine novae have been used to examine the properties of novae in carly-typc galaxies, 
measure the nova rate in M49, and assess the potential of novae as distance indicators. The 
main results of our study are as follows: 

• Compared to the M31 and Galactic samples, M49 may be under- abundant in slow, 
faint novae. Moreover, the distribution of novae decline rates in M49 is statistically 
indistinguishable from that observed for LMC Bearing in mind the small sample of 
novae on which our discussion is based, the M49 results seem to argue against a simple 
classification of novae in a bright, fast, disk population (which should be prevalent in 
the LMC) and a faint, slow, bulge population (to which all of the M49 novae should 
belong) . 

• At a distance modulus of 31.06 ± 0.10 mag, measured both using SBF and GCLF, the 
zero point of the Maximum Magnitude versus Rate of Decline relation for the M49 
novae is consistent with that derived from a sample of two dozen Galactic novae, with 
distances determined using expansion parallaxes. The agreement between the M49 
and M31 MMRD relations is less satisfactory, possibly owing to the large uncertainty 
associated with M31 internal extinction (which affects the maximum magnitude of the 
novae observed in this galaxy). In both cases, there seems to be a substantial difference 
in the shape of the MMRD relation in M49, the Milky Way and M31. 

• The mean magnitude of the M49 novae 15 days after maximum is marginally consistent 
only with one of three proposed calibrations based on Galactic novae. Furthermore, 
the magnitudes of the M49 novae seem to display a smaller scatter around maximum 
light than at 15 days past maximum. Altogether, these results caution against an 
indiscriminate use of novae as distance indicators. 

• The global nova rate in M49 is 77 = lOO^gg ycar"^, corresponding to a luminosity-specific 

nova rate vk in the range 1.7-2.5 year~^ IO^^^Lk© (depending on the adopted estimate 
for the X-band luminosity of the galaxy). This estimate accounts for observational 
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incompleteness, due both to the magnitude detection hmits, and to the selection criteria 
adopted in the detection of the variable stars. The value of measured for M49 is 
inconsistent with the predictions of the theoretical models, unless global differences 
are invoked between the novae progenitors in M49 and the Milky Way (against which 
the models arc calibrated). The luminosity specific nova rate in M49 is fully consistent 
with that measured in all other galaxies for which data are available, with the possible 
exception of the LMC. 

• Last but not least, the valuable lesson learned from our program is that, overall, ob- 
taining reliable light curves for novae is not a trivial task. Our program consumed 24 
orbits of HST time and lead to the discovery of nine novae. For comparison, 16 orbits 
of HSTimve were used to discover 52 Cepheid variables and measure an 8% distance to 
MlOO, also in Virgo (Ferrarese et al. 1996). In retrospect, a few changes to our observ- 
ing strategy would have been advisable. Color information would have been desirable 
at all epochs, and a two-day interval between subsequent exposures over the entire se- 
quence would have aided in the measurement of the novae light curve parameters. For 
new programs, the Advanced Camera for Surveys (ACS) would be more suitable than 
WFPC2 both because of the smaller pixel size (reducing the background noise) and 
higher sensitivity. Although these changes would lead to a better characterization of 
the novae light curves, they would entail a large program, likely requiring many dozens 
of HST orbits per galaxy. Even then, the low luminosity-specific nova rate, and the 
apparently large scatter in the MMRD and Buscombe-de Vaucouleurs relations would 
ultimately limit the usefulness of novae as distance indicators. SBF has been proven to 
be a rehable — and efficient — indicator for early type galaxies, while the GCLF has 
the potential of becoming one: they both seem more worthwhile choices for measuring 
distances. 
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Fig. 1. — Digitized sky survey image of M49 showing the location and orientation of our 
WFPC2 fields. The image measures 6' x 6'. 



Fig. 2.— HST/WFPC2 image of M49. This F555W image shows the WF2 field, with the 
positions of novae #1-5 indicated. 



Fig. 3.— HST/WFPC2 image of M49. This F555W image shows the WF3 field, with the 
positions of novae #6 and 7 indicated. 



Fig. 4.— HST/WFPC2 image of M49. This F555W image shows the WF4 field, with the 
positions of novae #8 and 9 indicated. 



Fig. 5. — The location of the nine novae discovered in M49, superimposed on an isophotal 
contour of the galaxy. The isophotes are drawn at 0.5 mag/arcsec^ intervals, from fiy = 16 
mag/arcscc^ to fiy = 20.5 mag/arcsec^. The WFPC2 footprint is also shown. The image is 
2".7 on the side, and has the same orientation as Figure 1 (N at the top and E to the left). 



Fig. 6. — Appearance of novae #1-5 at five different epochs. Each panel shows a 2'.'7 x 2'.'7 
region centered on the nova, in the F555W bandpass. Epoch numbers are indicated in each 
panel. 



Fig. 7. — Appearance of novae #6-9 at five different epochs. Each panel shows a 2'/7 x 2'/7 
region centered on the nova, in the F555W bandpass. Epoch numbers are indicated in each 
panel. 
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Fig. 8. — Light curves for novae #1, 2 and 3 in y (filled symbols) and / (open symbols). 
The lower and upper arrows show respective upper limits on the V and / magnitudes at 
epochs where the novae were undetected. 
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Fig. 9. — Same as Figure 8 except for novae #4, 5 and 6. 
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Fig. 10. — Same as Figure 8 except for novae #7, 8 and 9. 
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Fig. 11. — A comparison between the cumulative distribution of the novae detected in M49 
and the cumulative fraction of the underlying galaxy band light, normalized to the total 
light within the WFPC2 field. As is evident from Figures 2-5, the novae appear to be more 
centrally concentrated than the galaxy light. 
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Fig. 12. — The cumulative fraction of M31, LMC, Milky Way, and M49 novae as a function 
of decline rate, defined as = 2/^2, where t2 is the time it takes for the nova to decline by 
2 magnitudes (in F-band) after reaching maximum hght. The LMC and M49 samples are 
statistically indistiguishable. 
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Fig. 13. — Globular cluster luminosity function, 0(V^), for M49. Filled symbols indicate the 
luminosity function of Larsen et al. (2001). Open symbols show the luminosity function 
determined from our deep V and / images, scaled upwards by a factor (661/389) ~ 1.7. 
The solid curve shows an unconstrained Gaussian fit to the luminosity function derived from 
our observations; the apparent magnitude of the turnover is = 23.87 ± 0.06. The 

dashed curve shows the best-fit obtained when the apparent turnover magnitude is fixed at 
= 23.80; this is the value expected on the basis of the SBF distance modulus, (m — M)o 
= 31.06±0.10 (Tonry et al. 2001), and the absolute value of = -7.33 ± 0.04 given in 
Harris (2001). 
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Fig. 14. — The Maximum Magnitude versus Rate of Decline (MMRD) relation for the 
M49 novae, shown as circled symbols. Filled and open squares show M31 novae judged by 
Capaccioli et al. (1989) to have light curves of good and fair quality, respectively; their 
apparent magnitudes have been scaled to the distance and reddening of M49. The thick 
dashed curve shows the best fit MMRD relation for M31 (from Capaccioh et al. 1989, again 
shifted to the distance of M49), with la uncertainties shown by the thin dotted curves. 
Absolute magnitudes are shown at right under the assumption that (m — M)o = 31.06±0.10 
mag and A(y) — 0.07 mag for M49. The lower panel shows the deviations of the M49 novae 
from the M31 MMRD relation. 
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Fig. 15. — Same as Figure 14, except the comparison is now between the M49 (circled 
symbols) and Galactic novae (filled squares). The dashed and solid curves are the best 
fit "S-shaped" and power law fits to the Galactic novae from Downes & Duerbeck (2000); 
residuals between the former and the M49 novae are shown in the lower panel. 



-37- 




10 20 30 

t (days past maximum) 

Fig. 16. — (Upper Panel) Light curves for the five M49 novae having measured Knax and tmax- 
The hght curves have been shifted so that the times of maximum hght ahgn. Observed values 
are shown as filled symbols; open symbols indicate our best estimates for Knax- The novae 
have Vi5 = 24.77±0.19 at 15 days past maximum, corresponding to Mv,i5 = — 6.36±0.19. The 
standard deviation about this mean value is a ^ 0.43. (Lower Panel) The mean magnitude 
averaged over the five novae plotted in the upper panel, as a function of time following 
maximum. 
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Fig. 17. — Monte Carlo simulations of the global nova rate in M49. The histogram gives 
the number of times — based on 1000 simulations at each r) — that the simulated datasets 
contained exactly nine novae that would have been discovered using our observing strategy 
and data reduction procedures. The solid and dashed lines show our best estimates for the 
global nova rate and its l-cr uncertainties: r] = IOOI3Q year~^. 
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Fig. 18. — Luminosity-specific nova rates, vk-, for all galaxies with measured global nova 
rates. The global nova rates have been normalized by the total i^'-band luminosity of each 
galaxy, derived from 2MASS (upper panel - the open circles denote galaxies for which 2MASS 
magnitudes are not available, and i^-band magnitudes were derived from the RC3 5-band 
magnitudes as in the lower panel) or inferred from the RC3 5-band magnitudes via a color 
correction (lower panel - see text for further details). M49 is indicated by the circled point. 
For M33, two different — and inconsistent — estimates of the nova rate are available (see 
text for details). The dashed line shows the weighted mean, exluding M33: vk — 1.58 ± 
0.16 year~^ 1Q~^'^Lk,q in the lower panel, v'k — 2.41 ± 0.27 year~^ lO~i°Lx,0 in the upper 
panel. 
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Table 1. Log of Observations GO-8677 



Epoch Observation Date 



1 


2001 Apr 9 


1 


2001 Apr 9 


2 


2001 Apr 11 


2 


2001 Apr 11 


3 


2001 Apr 13 


3 


2001 Apr 13 


4 


2001 Apr 15 


4 


2001 Apr 15 


5 


2001 Apr 17 


5 


2001 Apr 17 


6 


2001 Apr 19 


6 


2001 Apr 19 


7 


2001 Apr 21 


7 


2001 Apr 21 


8 


2001 Apr 23 


8 


2001 Apr 23 


9 


2001 Apr 25 


9 


2001 Apr 25 


10 


2001 Apr 27 


10 


2001 Apr 27 


11 


2001 May 1 


11 


2001 May 1 


12 


2001 May 3 


12 


2001 May 3 


13 


2001 May 5 


13 


2001 May 5 


14 


2001 May 10 


14 


2001 May 10 


14 


2001 May 10 


14 


2001 May 10 



Julian Date^ Dataset Name 
(2,400,000+) 



52008 


.39844 


U6630101R 


52008 


.41406 


U6630102R 


52010 


.47266 


U6630201R 


52010 


.48438 


U6630202R 


52012 


.47656 


U6630301R 


52012 


.49219 


U6630302R 


52014 


.48438 


U6630401R 


52014 


50000 

• W w w w w 


U6630402R 


52016 


.49219 


U6630501R 


52016 


.50781 


U6630502R 


52018 


50000 


U6630601R 


52018, 


.51562 


U6630602R 


52020 

T_/ \J \J . 


.50391 


U6630701R 


52020 


.52344 


U6630702M 


52022 

\j \j . 


.44141 


U6630801R 

\J\JKJ\JyJ\J \j 


52022 


.46094 


U6630802R 

V-/ \J\JKJ\JyJ\J \j 




/I cqi O 


uooouyuiit 


52024 


.47266 


U6630902R 


52026 


.46094 


U6631001R 


52026 


.47656 


U6631002R 


52030 


.47656 


U6631201R 


52030 


.48828 


U6631202R 


52032 


.41016 


U6631301R 


52032 


.42578 


U6631302R 


52034 


.41797 


U6631401R 


52034 


.43359 


U6631402R 


52039 


.57422 


U6631501R 


52039 


.58594 


U6631502R 


52039 


.63281 


U6631503R 


52039 


.64844 


U6631504R 



Exposure Time Filter 



(seconds) 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1200 


F555W 


1300 


F814W 


1300 


F814W 
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Table 1 — Continued 



bpocn 


Observation Date 


Julian Date 
(2,400,000+) 


Dataset JName 


Exposure lime 
(seconds) 


b liter 


15 


2001 


May 


15 


52044 46094 


U6631601R 


1200 


F555W 


15 


2001 


May 


15 


52044.48047 


U6631602M 


1200 


F555W 


15 


2001 


May 


15 


52044.52734 


U6631603R 


1300 


F814W 


15 


2001 


May 


15 


52044.54297 


U6631604R 


1300 


F814W 


16 


2001 


May 


20 


52049.54297 


U6631701R 


1200 


F555W 


16 


2001 


May 


20 


52049.56250 


U6631702R 


1200 


F555W 


16 


2001 


May 


20 


52049.60938 


U6631703R 


1300 


F814W 


16 


2001 


May 


20 


52049.62500 


U6631704R 


1300 


F814W 


17 


2001 


May 


25 


52054.42578 


U6631801R 


1200 


F555W 


17 


2001 


May 


25 


52054.44531 


U6631802R 


1200 


F555W 


17 


2001 


May 


25 


52054.49219 


U6631803R 


1300 


F814W 


17 


2001 


May 


25 


52054.50781 


U6631804R 


1300 


F814W 


18 


2001 


May 


30 


52059.50781 


U6631901R 


1200 


F555W 


18 


2001 


May 


30 


52059.52344 


U6631902R 


1200 


F555W 


18 


2001 


May 


30 


52059.57812 


U6631903R 


1300 


F814W 


18 


2001 


May 


30 


52059.58984 


U6631904R 


1300 


F814W 


19 


2001 Jun 3 


52063.45312 


U6631101R 


1200 


F555W 


19 


2001 Jun 3 


52063.47266 


U6631102R 


1200 


F555W 



'^Modified Julian date at the middle of the exposure. 
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Tabic 2. Coordinates of Novae in M49 



Nova 


a(J2000) 


5(J2000) 


R 








(arcsec) 


1 


12:29:45.286 


7:59:42.73 


28.79 


2 


12:29:45.087 


7:59:39.95 


32.84 


3 


12:29:45.116 


7:59:54.95 


25.26 


4 


12:29:46.102 


7:59:46.61 


17.74 


5 


12:29:44.484 


7:59:45.93 


37.19 


6 


12:29:46.279 


8:00:37.51 


36.78 


7 


12:29:46.428 


8:00:15.15 


14.57 


8 


12:29:47.246 


8:00:20.35 


20.25 


9 


12:29:17.277 


8:00:23.18 


23.36 



Note. — Units of right ascension are hours, 
minutes, and seconds, and units of declination 
are degrees, arcminutes, and arcseconds. R is 
the distance of the nova from the center of 
M49. 



Table 3. Novae Photometry 



bpocn 


Julian Date 


Nova 1 




Nova 2 


Nova 3 


Nova 4 








V 


I 




V 


I 


V 


I 


V 


I 




1 


2452008 4075 


23.60±0.06 






22.64±0.04 




23.36±0.05 




>25.50 






2 


2452010 4818 


24.33±0.11 






23.13±0.04 




24.23±0.08 




>25.50 






3 


2452012 4887 


24.75±0.13 






23.67±0.08 




25.44±0.27 




>25.50 






4 


2452014 4964 


24.98±0.17 






23.90±0.06 




>25.90 




>25.50 






5 


2452016 5033 


25.74±0.26 






24.37±0.09 




>25.90 




>25.50 






6 


2452018 5102 


25.89±0.35 






24.91±0.21 




>25.90 




>25.50 






7 

8 


2452020 5179 
2452022 4575 


25.63±0.28 
25.39±0.24 






25.11±0.16 
25.28±0.17 




>25.90 
>25.90 




>25.50 
>25.50 




c 


9 


2452024 4644 


25.56±0.28 






25.61±0.19 




>25.90 




>25.50 






10 


2452026 4713 


>26.30 






25.69±0.25 




>25.90 




>25.50 






1 1 


2452030 4846 


25.78±0.30 






>26.30 




>25.90 




>25.50 






12 


2452032.4248 


25.74±0.33 






>26.30 




>25.90 




>25.50 






13 


2452034.4318 


25.37±0.21 






>26.30 




>25.90 




25.03±0.31 






14 


2452039.5863 


25.58±0.24 


>24 


80 


>26.30 


>24.80 


>25.90 


>24.60 


23.39±0.09 


22.95±0 


16 


15 


2452044.4706 


>26.30 


>24 


80 


>26.30 


24.43±0.26 


>25.90 


>24.60 


23.95±0.12 


23.07±0 


19 


16 


2452049.5547 


>26.30 


>24 


80 


>26.30 


>24.80 


>25.90 


>24.60 


24.14±0.15 


23.31±0 


21 


17 


2452054.4373 


>26.30 


>24 


80 


>26.30 


24.60±0.28 


>25.90 


>24.60 


24.28±0.13 


23.31±0 


21 


18 


2452059.5207 


>26.30 


>24 


80 


>26.30 


>24.80 


>25.90 


>24.60 


23.90±0.12 


23.39±0 


22 


19 


2452063.4651 


>26.30 






>26.30 




>25.90 




24.02±0.11 







Table 4. Novae Photometry 



Epoch 


T 1 • J. 

Juhan Date 


Nova 5 


Nova 6 


Nova 7 


Nova 8 






V 


/ 


V 


I 


V 


I 


V 


I 


1 


2452008 4075 


>27.00 




>26.30 




>25.20 




22.82±0.04 




2 


2452010 4818 


>27.00 




>26.30 




>25.20 




23.17±0.05 




3 


2452012 4887 


>27.00 




>26.30 




>25.20 




23.34±0.06 




4 


2452014 4964 


>27.00 




>26.30 




>25.20 




23.65±0.08 




5 


2452016 5033 


>27.00 




>26.30 




>25.20 




24.03±0.11 




6 


2452018 5102 


>27.00 




>26.30 




>25.20 




24.05±0.11 




7 

8 


2452020 5179 
2452022 4575 


>27.00 
>27.00 




>26.30 
25.82±0.83 




>25.20 
>25.20 




24.19±0.13 
24.39±0.13 




9 


2452024 4644 


>27.00 




23.59±0.10 




24.78±0.34 




24.45±0.14 




10 


2452026 4713 


>27.00 




24.41±0.15 




23.34±0.10 




24.90±0.23 




1 1 


2452030 4846 


>27.00 




24.27±0.10 




23.25±0.11 




>25.30 




12 


2452032.4248 


>27.00 




24.48±0.11 




23.58±0.10 




>25.30 




13 


2452034.4318 


>27.00 




24.57±0.12 




24.29±0.27 




>25.30 




14 


2452039.5863 


26.61±0.47 


>25.10 


24.87±0.15 


24.06±0.18 


24.77±0.34 


>24.00 


>25.30 


>24.60 


15 


2452044.4706 


23.91±0.04 


23.36±0.11 


24.95±0.17 


>24.80 


25.02±0.41 


>24.00 


>25.30 


>24.60 


16 


2452049.5547 


25.26±0.16 


24.45±0.23 


>26.30 


>24.80 


>25.20 


>24.00 


>25.30 


>24.60 


17 


2452054.4373 


25.38±0.17 


24.30±0.19 


>26.30 


>24.80 


>25.20 


>24.00 


>25.30 


>24.60 


18 


2452059.5207 


25.41±0.18 


24.37±0.20 


>26.30 


>24.80 


>25.20 


>24.00 


>25.30 


>24.60 


19 


2452063.4651 


25.38±0.19 




>26.30 




>25.20 




>25.30 
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Table 5. Novae Photometry 



Epoch Julian Date Nova 9 

V 7 



1 

J. 


9459008 4075 


>95 90 




2 


2452010.4818 


>25.90 




3 


2452012.4887 


>25.90 




4 


2452014.4964 


>25.90 




5 


2452016.5033 


>25.90 




6 


2452018.5102 


>25.90 




7 


2452020.5179 


>25.90 




8 


2452022.4575 


>25.90 




9 


2452024.4644 


>25.90 




10 


2452026.4713 


>25.90 




11 


2452030.4846 


23.99±0.12 




12 


2452032.4248 


23.77±0.08 




13 


2452034.4318 


23.82±0.08 




14 


2452039.5863 


24.48±0.14 


23.91±0.33 


15 


2452044.4706 


24.21±0.11 


23.90±0.35 


16 


2452049.5547 


24.96±0.21 


>24.50 


17 


2452054.4373 


25.11±0.21 


>24.50 


18 


2452059.5207 


25.18±0.24 


>24.50 


19 


2452063.4651 


25.70±0.40 
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Table 6. Properties of Light Curves for Novae in M49 



Nova 


^max 


Vmax 






t2 






(Julian Date) 


(mag) 




(mag) 


(days) 


(days^^) 


1 


<2452008.4075 


<23.60±0.06 


< 


-7.53 ±0.12 


7.78±0.62 


0.257±0.020 


2 


<2452008.4075 


<22.64±0.04 


< 


-8.49 ±0.11 


9.37±0.38 


0.213±0.009 


3 


<2452008.4075 


<23.36±0.05 


< 


-7.77 ±0.11 


4.48±0.40 


0.446±0.040 


4 


2452039.04±1.11 


23.33±0.13 




-7.80 ±0.16 


52.00±4.21 


0.038±0.003 


5 


2452041.42±1.95 


23.10±0.52 




-8.03 ±0.53 


7.53±1.14 


0.266±0.011 


6 


2452023.46±0.50 


23.33±0.28 




-7.80 ±0.30 


19.30±0.86 


0.104±0.005 


7 


2452027.33±0.50 


22.72±0.36 




-8.41 ± 0.37 


11.68±0.83 


0.171±0.012 


8 


<2452008.4075 


<22.82±0.04 


< 


-8.31 ±0.11 


17.71±0.81 


0.113±0.005 


9 


2452031.45±0.50 


23.74±0.22 




-7.39 ±0.24 


37.16±2.86 


0.054±0.004 



Table 7. Absolute and Normalized Nova Rates 



Galaxy T t^^ B K^mass {B - K)^'' A{By A{B)/ {m - Mf 

^K,B '^K,2MASS 







(year ^) 


(mag) 




(mag 


) 




(mag; 


) 


(mag) 


(mag) 




(mag) 




(year" 


-1 IO-ioLko) 


(year-i IQ-'^Lkq) 


LMC 


9 


2.5±0.5 


0.91±0.05 








2, 


.74±0. 


,10 


0.07 


0. 


,32 


18, 


.50±0. 


,13 


5 


.98±1.26 




SMC 


9 


0.3±0.2 


2.70±0.10 








2, 


.71±0. 


10 


0.24 


0. 


,13 


18, 


.99±0. 


,05 


2 


.50±1.68 




M33 


6 


4.6±0.9 

< 0.45 


6.27±0.03 


4, 


.11±0. 


,04 


2, 


.87±0. 


,10 


0.33 


0. 


,18 


24, 


.64±0. 


,09 


4 


.27±0.89 
< 0.41 


12.6±2.7 
< 1.22 


iVl i U i 





1 O 1 A 


o Q1 _i_n no 


0, 


.OiitU, 


DO 


Q 

O, 


,Z4±U. 


1 1 
i i 


U.UO 


0, 


,04 


oo 

zy. 


q /I 1 n 
,o4ltU. 




i 


.UiztU.oO 


i.OOltU.O ( 


M51 


4 


18±7 


8.96±0.06 


5, 


.05±0. 


,03 


3 


,43±0. 


10 


0.30 


0, 


,16 


29, 


,42±0. 


,27 


1, 


.53±0.61 


1.44±0.67 


MlOO 


4 


25±12.5 


10.05±0.08 


6, 


.59±0. 


,04 


3 


,84±0. 


20 


0.10 


0, 


,11 


31, 


,04±0. 


,09 


1 


.12±0.57 


l.|90±0.97 


M31 


3 


29±4 


4.36±0.02 


0, 


.98±0. 


,02 


3, 


,85±0. 


10 


0.67 


0. 


,35 


24, 


,42±0. 


,10 


1, 


.45±0.28 


5^5±0.88 
2.73±0.94 


M81 


2 


24±8 


7.89±0.03 


3 


.83±0. 


,02 


3 


.99±0. 


10 


0.34 


0. 


,35 


27, 


.80±0. 


,08 


1 


.64±0.56 


NGC5128 


-2 


28±7 


7.84±0.06 


3 


.94±0. 


,02 


3 


.38±0. 


11 


0.00 


0. 


,50 


28, 


.12±0. 


,14 


2 


.84±0.76 


2.67±0.75 


NGC1316 


-2 


135±45 


9.42±0.08 


5, 


.59±0. 


,02 


4, 


.15±0. 


,20 


0.00 


0. 


,09 


31 


.66±0. 


,17 


1 


.61±0.55 


2.33±0.86 


M87 


-4 


91±34 


9.59±0.04 


5, 


.81±0. 


,02 


4, 


.17±0. 


10 


0.00 


0. 


,10 


31 


.03±0. 


,16 


2 


.20±0.83 


3.43±1.38 


VirgoEs 


-4 


160±57 


9.46±0.10 








4, 


,26±0, 


11 


0.00 


0, 


,09 


31 


,17±0, 


,09 


2, 


.79±1.03 




M49 


-5 




9.37±0.06 


5, 


.40±0. 


,03 


4, 


.30±0. 


10 


0.00 


0. 


,09 


31, 


.06±0. 


,10 


1 


.71±0.61 


2.52±0.91 



''Nova Rate References: LMC: Capaccioli et al. (1990). SMC: Graham (1979). M33: Delia Valle et al. (1994) and Sliarov (1993). MlOl, 
M51: Shatter et al. (2000). MlOO: Ferrarese et al. (1996). M31: Capaccioli et al. (1989). M81: Moses & Shatter (1993). NGC5128: 
Ciardullo ct al. (1990). NGC1316: Delia Valle & Gilmozzi (2002). M87: Shatter et al. (2000). Virgo: Pritchet & van den Bergh (1987). 
M49: this paper. 

^Color References: LMC, SMC, M33, MlOl, M51, M31, M81, NGC5128, M87, Virgo: Shaftcr et al. (2000). MlOO: Aaronson (1978). 
NGC1316: Delia Valle (2002). M49: {B - V)o from NED and {V - K)q from Frogel ct al. (1978). 

^Errors on the internal and Galactic extinction estimates are assumed to be 20% and 16% respectively (the latter from Schlegel et al. 
1998). 

<^Distance References: LMC, SMC, M33, MlOl, M31, M81, MlOO: Ferrarese et al. (2000). M51, NGC5128, NGC1316, M87, M49, Virgo: 
Tonry et al. (2001). 



